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The radiation from an electron in a homogeneous magnetoplasma has some unusual properties 
as a consequence of the dispersive anisotropic nature of the medium. Attention is confined to emis- 
sion in the ordinary (whistler) mode frequency band below the cyclotron resonance and the extraor- 
dinary mode frequency band around the plasma frequency where the indices of refraction are appre- 
ciably greater than one and vary significantly. Due to the large indices, electrons can emit Cerenkov 
radiation over a limited band of nonrelativistic energies. The cyclotron radiation which is generated by 
the gyrations of electrons is complicated also by this property of the medium which permits both normal 
emission due to "slower than light" motion and anomalous emission due to "faster than light" motion. 
In the ordinary mode, for example, the anomalous cyclotron radiation is emitted into the forward hemi- 
sphere with respect to the guiding center motion of the electron whereas the normal radiation is emitted 
into the backward hemisphere. In this paper the frequency spectra and angular patterns of the average 
radiated power are calculated by the Hamiltonian method which avoids a direct calculation of the com- 
plicated electromagnetic field vectors. The theory of emission in dispersive anisotropic media with 
a hermitian dielectric tensor by Kolomenskii and Eidman is thoroughly reviewed and extended to 
include relativistic energies; the complicated analytic formulas for the power are evaluated for several 
special cases; and the results are applied to recent interpretations of very low-frequency (VLF) and 
low-frequency (LF) emissions from electrons in the magnetosphere. The main conclusions of the 
work are as follows: (1) In the ordinary (whistler) mode most of the energy is radiated along wave nor- 
mals at large angles to the magnetic field at frequencies other than the rectilinearly Doppler-shifted 
fundamental cyclotron harmonic which is contrary to assumptions of certain VLF emission theories. 
(2) The resonance singularities in the indices for a cold, collisionless plasma must be eliminated to 
achieve finite power levels, but unfortunately the dielectric tensor for thermal motion is extremely 
complex and for collisions is non-hermitian; consequently, an arbitrary upper limit is imposed on 
the indices in order to make a quantitative estimate of the power. (3) Based on this approximation 
the total power in the ordinary mode is a slowly varying function of frequency and electron energy with 
an average level of 10" 30 W/(c/s) per electron. (4) This level is inadequate to explain observed VLF 
signals on the basis of incoherent emission, but coherent emission from bunches of electrons can 
give the observed power level of 10 -14 W/cm 2 (c/s) above the ionosphere; hence, the onus of explaining 
the complex dispersion patterns of VLF emissions is left to the coherence mechanism. (5) The radi- 
ation in the extraordinary mode varies considerably with frequency and energy but an average power 
level is on the order of 10~ 25 W/(c/s) per electron which still requires some coherence to generate the 
observed level above the ionosphere; however, this energy cannot penetrate the ionosphere to account 
for the dispersion observed by ground-based receivers. 

1. Introduction 

The problem of radiation from an electron in a plasma in a homogeneous magnetic field occurs 
in astrophysics, radioastronomy, thermonuclear physics, and magnetospheric physics. Due to 
the dispersion and anisotropy of the magnetoplasma, the emission from an electron moving along 
a helical trajectory has some very unusual properties. The analysis is restricted to those fre- 
quency bands in the ordinary and extraordinary modes where the indices are appreciably greater 
than one. The polarization induced in the medium by the electronic charge produces Cerenkov 
radiation when the local phase velocity is less than the particle velocity, and acceleration due to 
the orbit gyrations of the charge produces cyclotron radiation which is appreciable in the low-order 
harmonics at Doppler-shifted frequencies. The cyclotron radiation consists of a "normal" emis- 
sion in which the source moves slower than the local phase velocity and an "anomalous" emission 
in which the source moves faster. 

The general nonrelativistic formulas for the frequency spectra and polar diagrams of the power 
radiated by gyrating electrons in a magnetoplasma have been given by Eidman [1958]. His analysis 
is based on earlier work by Kolomenskii [1953] in which the fundamental Hamiltonian method of 
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calculating radiation from a source is generalized for the case of an anisotropic dispersive medium 
with a Hermitian dielectric tensor. Although the formulas are available, the emission properties 
have not been analyzed appreciably because extensive numerical evaluation is required. 
P. S. Johnson [1962] used the expressions to evaluate the frequency spectrum of Cerenkov radiation 
from an electron moving parallel to the magnetic field. The equation relating the frequency and 
angle of emission has been solved numerically by Ellis [1962] for certain parameters appropriate 
for radiation from Jupiter. The power formulas were also evaluated by Bazhanova and Shafranov 
[1963] for the special case of normal emission in the ordinary mode from the fundamental cyclotron 
harmonic. However, the general problem of Cerenkov radiation and normal and anomalous cyclo- 
tron radiation from several harmonics has not been considered for electrons in helical orbits. 

Some approximate formulas have been derived independently. Several limiting cases of 
emission in a thermal magnetoplasma have been given by Pakhomov, Stepanov, and Aleksin [1962], 
and Pakhomov and Stepanov [1963a, 1963b]. Another method was employed recently by McKenzie 
[1963] to obtain the frequency spectra and polar diagram for Cerenkov radiation from a uniformly 
moving point charge. 

This paper has three objectives. First, the complete theory by Eidman [1958] and Kolomenskii 
[1953] is reviewed in order to include relativistic effects and pick up certain errors and omissions in 
the original analysis. Second, a numerical study of the formulas for Cerenkov and cyclotron radia- 
tion is presented to demonstrate the unusual properties of the emissions. Third, the results are 
applied to recent interpretations of the very low-frequency (VLF) and low-frequency (LF) emissions 
from the magnetosphere (3 to 300 kc/s). 

As a first approximation, the medium is represented by a cold, collisionless electron plasma. 
However, at the resonance frequency where the index of refraction becomes infinite, the power is 
also infinite. By introducing additional physical effects such as thermal motion or collisions, 
the singularity is removed but the expressions become so complicated that numerical evaluation 
is very difficult. As a compromise, an upper bound is imposed on the index which cuts off the 
emission at a finite level. Fortunately, in many cases this cutoff model is a good approximation 
to the index when thermal and collisional effects are small, and the physical properties of the power 
are described quantitatively in some parameter regions. 

The emission at high frequencies where the propagation is like that in a vacuum is not dis- 
cussed in this paper. Attention is confined to frequencies near the cyclotron and plasma frequen- 
cies of the electrons where the indices are most interesting. In the ordinary (whistler) mode the 
emission is confined to frequencies below the cyclotron frequency. In the extraordinary mode the 
emission is confined to a narrow band of frequencies near the plasma frequency. The index of 
refraction for these modes is generally greater than unity, and the electric vector of the wave ro- 
tates in the same sense as the electrons gyrate about the field (right-hand polarization). 

Due to the symmetry of the medium, the energy is radiated along cones about the magnetic 
field direction. The Cerenkov radiation is emitted entirely in the forward hemisphere with respect 
to the guiding center motion of the electron (velocity component along the static magnetic field). 
The normal cyclotron radiation is emitted entirely into the backward hemisphere in the ordinary 
mode, but into either or both hemispheres in the extraordinary mode. The anomalous cyclotron 
radiation is emitted entirely into the forward hemisphere. In general, the maximum power is 
radiated at an appreciable angle to the axis of symmetry. 

The anomalous cyclotron radiation is a consequence of the anomalous Doppler effect which 
was discussed some time ago by Franck [1942]. Since this property only arises when the source 
moves faster than the local phase velocity of the medium, the physical process is analogous to the 
Cerenkov mechanism; the charge induces a polarized wake in the medium which radiates electro- 
magnetic energy. An elementary example of the anomalous effect occurs in the acoustic case 
where the source moves faster than the speed of sound, and a stationary observer receives both 
increased and decreased (Doppler-shifted) frequencies when the source is receding but no signal 
when the source is approaching. 
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Finally, the feasibility of Cerenkov or cyclotron radiation from electrons as the source of VLF 
and LF emissions in the magnetosphere is discussed. A review of the early observations and the- 
ories and a classification system for the various emissions has been prepared by Gallet [1959]. The 
frequency-time patterns range from continuous wide bands of noise called hiss to discrete patterns 
with descriptive class titles: hooks, risers, quasi-vertical ("dawn chorus"), etc. The intensity of 
the wide band noise has been measured simultaneously at several frequencies by Dowden [1962a], 
and after ionospheric absorption effects were removed it was concluded that the source spectrum 
is relatively flat at a level of the order of 10~ 14 W/cm 2 (c/s). The morphology of VLF emission events 
has been recorded by Watts, Koch, and Gallet [1963] using an instrument which permits continuous 
observations. Recently, Helliwell [1963] has found that certain periodic VLF emissions arc trig- 
gered by whistler echoes. 

In addition to these ground-based observations, VLF emissions have been observed above the 
ionosphere by satellites. The VLF receiver aboard the Alouette satellite was operated by Barring- 
ton and Belrose [1963]. Their telemetered recordings of wide band noise exhibit a systematic 
decrease in the lower cutoff frequency as the satellite latitude increases. More detailed VLF 
emission data was obtained by Gurnett [1963] with instrumentation aboard the Injun III satellite. 
In addition to the telemetered wide band signal, a six channel frequency spectrum analyzer meas- 
ured the absolute amplitude of the VLF magnetic component. The rms amplitude at the satellite 
varied from the threshold of 10~ : * to a maximum of 5 X 10 -2 gamma ( ~ 10 -14 W/cm 2 (c/s)). On 
several occasions complex VLF emissions were triggered by spheric impulses. A clear correla- 
tion among VLF emissions, auroral enhancements, and electron precipitation was also revealed. 
In a rocket flight to 1700 km, Walsh, Haddock, and Schulte [1963] observed intense radio noise 
near 1 Mc/s, which they attribute to emissions in both modes. 

Several mechanisms have been suggested as the source of VLF emissions. Gallet and Helli- 
well [1959] proposed that geomagnetic field alined streams of electrons interact, as in the traveling 
wave tube mechanism, to produce amplification of existing background electromagnetic radiation. 
This mechanism was recently considered in greater detail by Dowden [1962b] who showed that 
traveling wave tube amplification of selected bands of VLF noise can explain several character- 
istics of hiss. The possibility of Doppler-shifted VLF cyclotron radiation from protons was pro- 
posed by MacArthur [1959] and later extended by Murcray and Pope [1960a, 1960b]. Dowden 
[1962c, 1963] also showed that Doppler-shifted cyclotron emissions from discrete bunches of 
electrons produce hook dispersion patterns which fit the experimental data quite well. Bell and 
Buneman [1964] found a plasma instability in the interaction between a whistler-mode wave and 
a gyrating-electron stream which may initiate coherent triggered emissions. Triggering mechan- 
isms and cyclotron emissions are also discussed by Hansen [1963]. The transfer of energy 
between whistler-mode signals and energetic charged particles has been examined for several 
VLF emission mechanisms by Brice [1964]. 

In most of these theoretical analyses only the one-dimensional longitudinal equations (parallel 
to the geomagnetic field) were considered so that wave normal directions in which the maximum 
Cerenkov and cyclotron power is available were tacitly omitted. Furthermore, no estimates of 
the total power or its spectrum were attempted. With the power estimates presented here it is 
demonstrated that incoherent emission is undoubtedly inadequate, but coherent emission from 
bunches of electrons may provide the observed power level. 

2. Hamiltonian Method 

The Hamiltonian method for calculating the energy radiated by a charged particle in an ani- 
sotropic dispersive medium is described in some detail. The solution which was formulated by 
Kolomenskii [1953] is a generalization of a method employed by Heitler [1954 and earlier editions] 
and others for calculating particle emissions in free space. The name of the method arises because 
the coefficients in the eigenwave expansion for the vector potential of the radiation field are just 
the canonical variables of the Hamiltonian for the system. 
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Any homogeneous medium which is electromagnetically dispersive and anisotropic is char- 
acterized by a complex dielectric tensor ^(co, 0, $), where o> is the angular propagation frequency 
and 6 and (f) are the polar and azimuthal angles measured from some reference direction. It is 
necessary for the following analysis to assume that this tensor is Hermitian 

4=e*, (2.1) 

where the asterisk denotes the complex conjugate and transpose. Furthermore, the magnetic 
susceptibility of the medium is assumed to vanish so that the permeability /ul of the medium is that 
of free space. The implications of these assumptions in a magneto-ionic medium are discussed 
in the next section. 

A charged particle moving through such a medium may radiate electromagnetic energy if it 
is being accelerated or if its velocity is greater than the phase velocity of the medium. The 
electric and magnetic fields which are generated by a particle with charge e, trajectory r e and 
velocity v are described by Maxwell's field equations (Gaussian units) and some familiar subsidiary 
relations: 

VXH-Dlc = 47r(e/c) v 8(r - r e ) (2.2) 

VxE+B/c-0 (2.3) 

V-B = (2.4) 

V • D = 47re8(r — r e ) n 5) 

D = e E,B = H (2.6) 

where 8 is the Dirac delta function and the dot denotes d/dt. 

By introducing the electromagnetic potentials, A and <I>, which are defined by the field 
quantitites, 

E=-V*-A/c, (2.7) 

B = VXA, (28) 

equations (2.2) to (2.6) reduce to 

V X (V X A) +€_• VO/c +£• A/c 2 = 4tt- (elc) v8(r- r c ), (2.9) 

- V • e_- V* — V • €_• A/c = 477-e8(r — r e ). (2.10) 

Since the potentials are not uniquely defined by (2.7) and (2.8) a Coulomb gage condition 
may be imposed on them: 

V-e/A^O (2.H) 

In this gage, 4> is just the static potential of the source charge and D si = €_• V4> is the corresponding 
static vector field. It will be shown that the electromagnetic radiation fields E em and B em are 
derivable from the vector potential A alone. 

In order to express these equations in canonical form it is necessary to consider only that part 
of the field which lies inside an arbitrarily large but finite cubic volume t=L 3 , where L > charac- 
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teristic wavelength. In this volume the potentials are represented by a superposition of eigen- 
functions A\ and <£\, 

A = 2gxWA x (r), (2.12) 

$ = 2gox(0^\(r), 

x 

where q\ and gox are expansion coefficients. The eigenfunctions are required to satisfy the 
boundary conditions, 

A A , VxAx, <J> A , Vd>x (2.13) 

periodic on the surface o~ of the volume r, 

and the homogeneous wave equations, 

V X (V X Ax) - cox 2 | • Axle' 1 = 0, (2.14) 

V 2 (V(D X ) +coo 2 xf • V4>x/c 2 = 0, (2.15) 

where the eigenvalues co\ and a> x are real constants. Since the divergence of a curl always van- 
ishes, the eigenfunctions Ax defined by (2.14) also satisfy (2.11) as required. The orthogonality 
of the set of Ax with respect to the weighting function e can be verified using (2.1), (2.13), and 
(2.14) and the divergence theorem, 



^-^) jdrA*'€'A x = c 2 fdr[Ax- V x V x A*- A* • V x V x A x ] 

= c 1 \ rfrV • [A* X V X Ax - A* X V X A M ] = 0, (X # fi). 



The orthogonality of the set of V^x with respect to e is verified analogously. The normalization 
of the eigenfunctions is chosen such that 



/ 
/ 



<irA* • c • Ax = 47rc 2 8x M , 

dr V<f>* • € • VOx = 4ttc 2 o Am . (2. 1 7) 



The orthogonality of Ax and V^> M for all X, fi with respect to the weighting function e can be 
proved with the aid of a vector identity, 

/ dr[h • V X (V X a) - V X b • V X a] =- fd<T- (b X V X a). (2.18) 

Inserting a = Ax and b — V^*, the surface integral vanishes due to the boundary conditions and 
the last term of the volume integral obviously vanishes to that with the aid of (2.1) and (2.14) 

JdrA* • f • V<JV= [/ drVO* • e • A x ]* =0. (2.19) 

The expansions for A and cj> permit (2.9) and (2.10) to be reduced to a system of uncoupled 
elementary scalar equations for </x and gox- If (2.9) is multiplied by A* and integrated over r, the 
orthogonality relations derived above reduce the vector equation to simple harmonic oscillator 
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equations for the set of q^ 

q„^^lq^={elc)\ • A*(r e ). (2.20) 

The forcing function depends on the dielectric tensor since A* must satisfy (2.14), and on the time 
through r e and v. Similarly multiplying (2.9) and (2.10) by V<£>* and integrating over r yields 

*=(e/c)vV$J(r e ), (2.21') 

c 2 ?0M = eO*(iv), (2.21) 

respectively. Since r e is explicitly time dependent, d<£>*{r e )ldt = \ • VO*(r P ), and therefore (2. 21') 
is redundant because it is just the time derivative of (2.21). Hence, the complete solution for the 
static and electromagnetic fields can be obtained in principle by solving (2.20) and (2.21) with the 
appropriate set of A M and <I> M for the medium. 

Since A and <£ are determined by uncoupled systems of equations, the electric and magnetic 
fields which they prescribe are also uncoupled. In particular, the electromagnetic fields are 
obtained from the vector potential alone using (2.7) and (2.8), 

£em=-A/c, B e m = VXA. (2.22) 

In the present application only the energy of the radiation field is required, and it is given by 
the general expression 

W = -tr \ dt ' J dT ^ E em ' Dem + B* m ' Hem). (2.23) 

The magnetic energy term is simplified with the aid of (2.18). Inserting a = A\ and b = A*, noting 
that the surface integral vanishes by the boundary conditions, and using (2.14), it is found that 

[rfrVxA*-Vx^ x =(co 2 x /c 2 )[^rA*- c • Ax. (2.24) 

Hence, using the normalization of (2.17), (2.23) reduces to 

W(t) = \df 2 (q*q + o>l€<lJ =ll (<?*<7m + «<7m) =2«V (2.25) 

J H. /i- M 

This expression is identical to that for the Hamiltonian of a system of harmonic oscillators which 
describe the radiation field. 

The polarization of the electromagnetic waves can take two arbitrary directions which corre- 
spond to the two possible modes of propagation in the medium. In the foregoing equations these 
different modes of polarization are denoted by different indices, X(or /x). In the following equations 
only one mode is considered. 

In order to evaluate the summation in (2.25) it is necessary to assume a particular form for the 
eigenvectors A M which satisfy (2.13), (2.14), and (2.17). Here, for convenience,- they are chosen 
to be 

A„ = V47T c(aJn^L 312 ) exp ( — ik M • r), (2.26) 

where a M is the polarization vector, n^ is the index of refraction, and k M is the propagation vector. 
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(The sign convention in the exponent is opposite to that used by Kolomenskii and Eidman, but it 
is essential for the equations to be consistent with the polarization convention used in magneto- 
ionic theory.) The orthogonality of A M in r is assured if the components of k M are given by k^-q 
= 2777WL where the fir) are integers and r\ —x, y, z. The components of a M and n M are determined, 
in principle, by (2.14) and (2.17) which reduce to 



(c 2 /^ 2 ) [A:£a M — k^(k M • a fJL )] = €_- a M , 



(2.27) 



Since the wavelength is small compared to L (L — > 0° ), /c-space is approximately continuous 
and jit-space can be represented by continuous spherical coordinates jh, 6, cp without loss of gen- 
erality. (The polar and azimuthal angles, 6 and (/>, are independent of the space variables in r.) 
Furthermore, by definition, conic = k = IttjiIL so that the frequency a> will be used as the inde- 
pendent variable with k(o), 0, (/>), n{co, 0, $), and a(a>, 0, 0), in general, and the subscript /x is omitted 
hereafter. 

The number of oscillators in r which radiate into the solid angle dfl = sin Odddqb and the fre- 
quency interval between co and to + dco for a prescribed polarization is given by 



rd^yu — jjL 2 djJidflT = 



uWL 3 

(2iTcf 



-, co dn\ 

1 -r - — I dcodilr = pdcodilT, 



(2.29) 



where p(w, 0, (f>) is the density function of the radiation. (In Kolomenskii's [1953] and Eidman's 
[1958] papers the frequency dependence of n was omitted in the definition of p.) 

With a continuous variable for /x the discrete summation in (2.25) transforms into a volume 
integral over the continuous variable w. Hence, the differential energy of the radiation field can 
be expressed as 



d 3 W = wpdo)dn, 



(2.30) 



where w (co, 6, c/>, t) is the energy of a forced oxcillator in dcodfl. 



3. Emission in a Magneto-Ionic Medium 

The foregoing general equations are developed further in this section for emission in a magneto- 
ionic medium which consists of a uniform static magnetic field and a homogeneous electrically 
neutral plasma. Such a medium is dispersive and anisotropic with respect to the direction of the 
magnetic field. The analysis follows Eidman's method with some modifications and corrections. 

The xyz coordinate system is oriented such that the z-axis is parallel to the static field B . If 
collisions are neglected, the orbit of a relativistic electron in the static field consists of a gyration 
about the field direction and a drift along it as shown in figure 1. Hence, the position (r e ) and 




Figure 1. The position and velocity coordinates of an electron on a 
helical trajectory. 
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FIGURE 2. The components of the electric vector E of the radiation 
field in a magnetoplasma. 

The wave normal k and the Poynting vector S are not parallel in general. 



velocity (v) coordinates along the x, y, z axes are, respectively, 

yn cos y^cont, yn sin y~ x oj B t, v 2 t, (3.1) 

— v\ sin y~ l o)Bt, V\ cos y _1 cofl£, f2, (3.2) 

where the subscripts 1 and 2 refer to components perpendicular and parallel to B and the rela- 
tivistic corrections have been introduced explicitly. The nonrelativistic cyclotron or gyrofrequency 
is related to B by a) B = |e|B /Mc, and the nonrelativistic gyroradius is given by ri = Vi/co B . The 
relativistic factor y is defined as usual by 

y = (l-p)-H\ p 2 = (v\ + vl)lc 2 . (3.3) 

The propagation vector k is at some angle 6 with respect to B , in general, and it is arbitrarily 
placed in the yz plane since there is azimuthal symmetry, 

k = kk = k sin 6'y + k cos 6 z, (3.4) 

where k y j» z are unit vectors. The components of the electromagnetic fields are taken along the 
orthogonal directions x, 0, k as shown in figure 2. 

The dielectric tensor for a magneto-ionic medium can have several forms depending on which 
properties of the plasma are important. In addition to the frequency domain, the effect of thermal 
motion, collisions, or ionic motion alter the dielectric tensor appreciably under certain conditions. 
Elementary derivations of the tensor for a cold plasma are given by Ratcliffe [1959], Ginzburg [1961], 
and others. Several papers have treated electromagnetic wave propagation in hot plasmas re- 
cently, and an extensive review is provided by Stix [1962]. It is noteworthy that the tensor is 
hermitian when thermal motion or ionic motion is included, but is not hermitian when collisions 
are important. 

The polarization of the electric vector of the wave field and the index of refraction for the me- 
dium are obtained from the components of the dielectric tensor. In terms of the spatial configura- 
tion in figure 2 the polarization coefficients are defined by 

E d /E x = iae((D, 0), E k /E^ = ia k {co, 0), (3.5) 

and the index is denoted by n(u>, 0), where a e and a k may be complex but n must be real in the 
following solution. 

The components of a (and n) in this medium are determined by (2.27) and (2.28) for a prescribed 
e. Since the electric field of the propagating wave must satisfy the wave equation, its amplitude 
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vector E must satisfy (2.27). Hence, the solution of these four complicated algebraic equations 
is obtained indirectly by assuming that a = #cE/£j?. The constant of proportionality is easily shown 
to be 

k = (1 + \E e IEA 2 )- il2 = (1 + a$)-«\ (3 6) 

(Eidman's normalization factor k= 1/ v2 does not satisfy the equations unless cos 0=±1.) There- 
fore, the components of a in the Cartesian coordinate system of the electron are given by 

a x = K, a y = iKa y , a z = iKa z , (3.7) 

where 

(Xy^ote cos + ak sin 0, 

(3.8) 

a z = ak cos — as sin 6. 

With the foregoing explicit expressions for A, i>, v, k, and a given by (2.26), (3.1), (3.2), (3.4), 
and (3.7), respectively, the forcing function in (2.20) takes the form 



(e/c)\ - A*= • v jm (—iv\ siny- l (x)Bt + a y V\ cos y~ 1 aj f it J r a z v 2 ) 

exp (+ ik sin yn sin y^wat -h ik cos 6 vtt). (3.9) 
This expression can be simplified by introducing the Bessel function expansion 

exp (+ iX sin 0) = ]T J 8 (X) exp (+ is </>), (3.10) 

— x 

so that 

[a y cos <\) — i sin (j>) exp (-f iX sin 4>) 

= ^2 (<*vJm + ctyJs-i+Js+i-Js-i) exp (+ u (/>) = 2 (-./.' + <W*W exp (+ is </>), (3.11) 

where the prime denotes differentiation with respect to X. Hence, (2.20) for the forced harmonic 
oscillator has the simple form 

+ x 

« q + Q) 2 q = ^C s exY>(+ip s t), (3.12) 

— x 

where 

p s (co, d) = s(t)B/y + V2k cos = scos/y + /Sa o>n(a>, 0) cos (3.13) 



and 

^4tt iecK 



Qa),«=- 



nL 3 / 2 



- j;(X) + (a^/3,/* + a,/8 2 )7,w] , (3.14) 
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with 

X(co, 6) = yr\k sin 6 = y/3ira(co, 0) (oj/cob) sin 6 (3.15) 

and 

Pl=Vi/c, p2 = V2lc. (3J6) 

With the initial conditions q = q = at t = 0, the general solution of (3.12) is 

£^ £2* Texp (ip s t) — exp (+ico*) . exp (ip s t) — exp ( — tco/)~| 
« = 2> = 2 C 4 2o,(a,-p.) + 2a,(a, + p.) J' < 3 - 17 ) 

The sum of the energies of all the modes then vanishes at £ = and is given by 

foo 1 



w = ^ y(qfqs + (D 2 qfq s ) + cross product terms 



— oo 



(3.18) 



^1 ,_ , 2 |" 1 — cos (w-p s )T 1 — cos (qj + p s )r i i 
^2 |Cs| L («-p s ) 2 (o)+ Ps ) 2 J ] 



— cos (a)—p s )T , 1 —cos (co + p s )T ~\ ^oscillating cross 

product terms 



at t = T. In the limit 7 — ■> oo, the explicit summation terms reduce to Dirac delta functions, 

T . 1 — cos (co± p s ) T sin (co±p s ) T , , 

Lim T~^ — \2T = Lim — w X \ =8(o)±.p 8 ),. (3 19) 

The first equality can be demonstrated by partial integration, and the second is discussed by 
Schiff [1955]. The average (unweighted differential) energy per unit time radiated into r by the 
electron is given by 

Lim w/r=| § |C S | 2 [8(w-p.) +8(ft>+p,)], (3.20) 

—oo 

where the oscillating transient terms vanish as 1/T. 

The limiting process which introduces the 8-functions is justified provided that the radiation 
never reaches the boundary <j of the volume r. This is assured because L and T can approach 
infinity and simultaneously satisfy the inequality L > cT/n. 

The mathematical interpretation of the 8-function terms depends on the propagation properties 
of the vector-potential quantity 

<7 M A M = q^ 0, cpMAo,, e, 4>(r) ~ exp i(±cot — k • r), (3.21) 

where q is given by (3.17), A is given by (2.26), and a) is always positive. If a particular set of 
values, wi, 5i ^ 0, d\, satisfy io = + p s , then by (3.13), coi, — Si, #1 + 77 satisfy a) =—p 8 . Therefore, 
the term 8(w— p s ) which is associated with the first term in (3.17) has a propagation exponent of 
the form -\-co\t — k\j sin 6\ — k\z cos 6i, where k = conic, and the term d(co-\-p s ) which is associated 
with the second term in (3.17) has a propagation exponent of the form — to\t + k\j sin 6\-\-k\Z 
cos 0i. Hence, these two terms describe identical radiation which is called "normal" emission 
because of its dispersion relation. Similarly, if CO2, 52 < 0, 62 satisfy o) = -\-p s , then 0)2, — 52, 02+tt 
satisfy co = — p s , and the two terms have a common propagation exponent of the form ±(co2t — k 2 y 
sin Qi — kzz cos 62). This radiation is called "anomalous" because of its unconventional dispersion 
relation which is discussed below. 
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Due to this duplication of emission properties, the energy expression can be simplified. By 
introducing the explicit property ao. a(cd, 6) = — ote, k((*>, O + tt), which is valid in the magneto-ionic 
medium, and using the properties of Bessel functions it is readily verified that C s (oj, 0) 
= — C-s((0, 0-\-tt). Therefore, the energies given by 8(oj — p s ) and 8{oj + p s ) are identical, and the 
radiation emitted by the electron is described completely by the terms 

Lim WsIT=tt\C s \ 2 8{(d — Ps), (3.22) 

T—><x> 

where 5 = 0, ± 1, ±2, . . . . The frequency 00 and polar angle 6 of the emission are related by 
the anisotropic dispersion equation a> = + p s , henceforth called the emission equation, which can 
be expressed in the form 

fe cos 8 = ((o — sa)Bly)l<on((o, S ) . (3.23) 

The physical explanation of the emission terms is provided by the properties of the emission 
equation. The s = term describes Cerenkov radiation, emission from induced polarization in 
the medium, because O must satisfy the usual relation cos Oo = l/fen. Since n is independent of 
the sign on the static-field vector-component, the Cerenkov emission is in the forward hemisphere 
with respect to the velocity component fe(fe cos O > 0). For s 7^ the emission is due to the 
cyclotron harmonics |s|cu#/y. The distinction between the normal and anomalous emission is 
obtained from the emission equation relation fen cos 6 S ^ 1 for s ^ 0. Hence, the guiding center 
moves slower than the parallel component of the phase velocity when s > and a normal Doppler 
effect is obtained, but it moves faster when 5 < and an anomalous effect is obtained. The 
normal emission appears in either the forward or backward hemisphere {fe cos X ^ 0) depending 
on whether co ^ soWy, whereas all the anomalous emission is in the forward hemisphere (fe cos 6 S 
> 0). Hence, frequencies below the source frequency |s|aWy are emitted into both hemispheres 
when the propagation modes permit. The simultaneous emission of both the normal and anoma- 
lous radiation is due to the dispersive anisotropic nature of the medium which allows waves of a 
given frequency to be reinforced in many directions at once. The general solution of (3.23) for 
cos 0s(<o) and its application to VLF and LF emissions in the magnetosphere is considered in the 
following section. 

The Cerenkov and cyclotron power which is radiated by the electron into dojdfl is given by 
terms of the form 

<PP S = Lim (wslT)pda)(m, 5=0, ±1, ±2, . . . (3.24) 

where the definition of average power and (2.30) have been used. In order to perform the inte- 
gration over frequency or solid angle, the delta function must be reduced to an elementary form 
by the expansion, 

8(f(x)) = 8[(x-xs)f(x s )] = [f'(x s )\-*8(x-x s \ (3.25) 

where x = or go and x 8 is that value of 6 or co for which o>=p s . Therefore, with p and C s given by 
(2.29) and (3.14), the power radiated into the linear frequency interval df=d(ol27T is obtained by a 
trivial integration, 

dPsldf= (27re*lcf3.>)a>{K> \rfeJ' 8 {X) + {ctysfelX + a z fe)J s (X)] 2 

• [l-\-((t)/n)dnld(o]l\l — (dnld6)ln tan 0\}e=6s(a>h (3.26) 
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where 6 s (oj) must satisfy (3.23). Similarly, the power radiated into the solid angle dfi is given by 
(PPsldil = (eV27Tc){(o 2 nK 2 [- pj',(X) + (otysfiilX + a z p 2 )MX)] 2 

• [l-\-(c*ln)dnlda)]l\l- p 2 n cos © [l + (a>/7i)dn/da>]|}«-«^), (3.27) 

where a), s (0) must satisfy (3.23). (These expressions for the power differ from those given by Eid- 
man [1958] in the factors p, k, and the sign on J' 8 , as well as including relativistic energies.) 

Two limiting cases of these power expressions are of special interest. First, consider /3i = 
which gives 

dPoldf=(27Te 2 /c/3 2 )oj{K 2 a 2 z l3 2 2 • [l + (a)ln)dnlda>]l\l-(dnldO)ln tan 0\}o=e<t») (3.28) 

for the Cerenkov radiation and the cyclotron terms vanish. Secondly, consider p 2 = 0, which 
gives ps = sa)Bly and 

(PPsldSl = (e 2 l27Tc){a> 2 nK 2 [- pj' s {X) + a y spj s (X)/X] 2 ■ [1 + (co/n)arc/dw] }„ =W y (3.29) 

for 5^1. Here only the normal cyclotron emission is present since the motion needed for the 
anomalous effect is absent. (Although the index is omitted from the emission equation in this 
case, the frequencies are restricted to the propagation bands of the modes, of course.) 

In free space where 71 = 1, a$=±l, afc = 0, the cyclotron terms in (3.29) reduce to the Schott 
formula [see Landau and Lifshitz, 9-8, 1951] when the power from both modes is added (5 > 
only since n= 1), 

d 2 Psld^ = {e 2 l2rrc){s 2 oy 2 ly 2 )\fiy' 2 {^s sin 0) + cot 2 6J*(p lS sin 0)] (3.30) 

where X from (3.15), k from (3.6), and co = so)Bly have been introduced explicitly. Since n > 1, 
a$ — 1, and ecu ^ 1 (for sin 6 > 0) in a magnetoplasma, a comparison of (3.30) and (3.29) shows that 
the magnitude of the power radiated into a plasma is far greater due to the interaction with the 
medium, but the emission equation limits the angles of emission (or frequency bandwidth). 

Although the power expression d 2 P s ldCl is evaluated at 6 and (j> (</> is arbitrary by azimuthal 
symmetry), this power follows a ray path along 0' and where 6' 7^ 6, in general. It has been 
shown by Scott [1950] that the ray path is in the direction of the time-averaged Poynting vector 
< S > = (c/47r) < E em X Hem >. With the electromagnetic fields denned by (2.22) and a propor- 
tional to E from (2.26), the ray path direction is along 

< S > a < E X (k X E) >. (3.31) 

Since E can have an appreciable component along k in an anisotropic medium, the electromagnetic 
energy does not propagate in the k direction, in general. For the specific as and otk given in the 
next section, it is readily verified that the direction of < S > is within the acute angle between 
k and ± B . 

4. Emission in a Cold, Collisionless, Electron Plasma 

In order to evaluate the allowed frequencies and radiated power from a relativistic electron, 
it is necessary to choose a model for the plasma. In this section the plasma is assumed to be so 
cold that thermal motion can be neglected, so tenuous that collisions can be neglected, and the 
frequencies of interest are sufficiently high that ionic motion can be neglected. Such a model is 
adequate for many physical systems including several aspects of VLF propagation in the magneto- 
sphere. In the present application, however, this elementary model is found to be unsatisfactory 
at certain frequencies and angles of emission. 
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Due to the complexity of the formulas and the large number of independent parameters, the 
behavior of the emission is demonstrated by several examples. In particular, the specific pa- 
rameter values are appropriate for trapped-electron orbits at sample points in the geomagnetic 
field. A brief summary of orbit and medium characteristics for the magnetosphere is given in 
the appendix. 

Since the propagation frequency is an independent variable here, it is convenient to introduce 
the dimensionless quantities F and P which are defined to be the propagation-cyclotron frequency 
ratio and the plasma-cyclotron frequency ratio, respectively, 



F = 0)/ (Ob , P = 0)pI0)b- 



(4.1) 



The plasma frequency is related to the electron density No by co 2 — 4?7rNoe 2 /m. In terms of the 
spatial configuration in figure 2 and these new variables, the dielectric tensor components for the 
medium described above have the form [Ginzburg, 1961] 

e^=l + P 2 /(l-F 2 ), 

e^=l-i >2 (sin 2 $-F 2 )IF 2 (l-F 2 ), 

ekk=l-P 2 (cos 2 0-F 2 )IF 2 (] -F 2 ), 

(4.2) 

€ xe = - e ex = + iP 2 cos 0/F\ 1 - F 2 ), 

e.rk = -e k . r = +iP 2 sin <9/F(l-F 2 ), 

e«i: = + e k „ = + P 2 cos 6 sin 0/F 2 (] -F 2 ). 
For these specific components, the indices of refraction for the medium are found to be 

n 2 ± = l + 2P 2 (P 2 -F 2 )/D ± , (4.3) 

with 

D ± =± [F 4 sin 4 6 + 4F 2 (P 2 -F 2 ) 2 cos 2 0]^ 2 -2F 2 (P 2 - F 2 )- F 2 sin 2 0, (4.4) 

and the polarization coefficients defined in (3.5) are given by 

a e ±=-F cos 0/ [P + /»/(/£ -1)] 

= - 2F(P 2 - F 2 ) cos 01 [- F 2 sin 2 6 ± (F 4 sin 4 6 + 4F 2 (P 2 ~ F 2 ) 2 cos 2 O) 1 ' 2 ], (4>5) 

a k± =-(n 2 ± -l)F sin 6/(P 2 - F 2 ) 

= - [P 2 F sin e-a e± P 2 cos sin 0]/[P 2 (cos 2 - F 2 ) - F 2 (l - F 2 )] . (4 6) 

In this notation the upper (+) sign denotes the ordinary (whistler) mode and the lower (— ) sign 
denotes the extraordinary mode of polarization (Ginzburg's nomenclature). The expressions for 
ae± and a k ± describe a wave of the form exp i{o)t — k • r). 

The indices given by (4.3) and (4.4) have several unusual properties. They have resonances 

(n 2 ± =±*>) at 



-(P 2 +l)-h-LV* --r ±j «*i— cus- uy- 



FL = -(P 2 +l)=h-[(P 2 -f 1) 2 -4P 2 cos 2 6] 1 ' 2 
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or 



cos 2 0fl = F 2 [l + (l-F 2 )/P 2 ] 



(4.8) 



as well as at F = 0, and they have cutoffs at 



1 1 J/ 2 

F = P (n + = 0) 9 F x , z = ±± + (± + P 2 ) (n- = 0). 



(4.9) 



From these relations it can be verified that F R+ < P and |cos 0|, F x > Fr- > P and 1, and F z < P. 
In a slowly varying medium, electromagnetic energy is absorbed at a resonance and reflected at a 
cutoff, in general. 

In the frequency band < F < Fr + , the ordinary mode index n + is real and appreciably greater 
than unity; for Fr+ < F < F , n+ is imaginary (evanescent); and for F < F < oo, n+ asymptotically 
approaches unity. In the frequency band < F < F z , the extraordinary mode index n~ is imag- 
inary; for Fz< F < F R -, ji- is real and increasing; for Fr- < F < F x , n~ is again imaginary; and 
for F x < F < oo, n- asymptotically approaches unity. 

At very high frequencies where the index is near unity, the radiation from an electron is 
described approximately by the free-space formulas derived by Schwinger [1949]. For those 
frequencies where the index is imaginary, emission of electromagnetic radiation cannot occur. 
Hence, attention will be confined hereafter to the ordinary-mode frequency band < F < F R +, 
where n+ > 1, and to the extraordinary-mode frequency band F z < F < F R -, where n- > 0. (Since 
ionic motion has been neglected, the frequencies are further restricted here to F >> mjM, Pm/M, 
where M is the ion mass.) 

The ordinary mode frequency band < F < Fr + is often called the whistler mode because VLF 
radio whistlers propagate in this mode. The whistler-mode index of refraction is plotted in figure 
3 as a function of F with parametric and P = 2. For a given F (and P) it is apparent that |cos 0\ 
> cos Or. This mode is frequently represented by its quasi-longitudinal approximation 



n%= l+P*IF(\co8 0\ -F) for F 2 sin 4 0/cos 2 « (P 2 -F 2 ) 2 , 



(4.10) 



and its resonance value, | cos 0|, is also indicated in figure 3. When the solution of the emission 
equation is near the resonance, the discrepancy between Fr+ and | cos 6 \ can cause serious errors 
in the power evaluation if the approximation is used. 
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Figure 3. The ordinary (whistler) mode index of refraction. 

The quasi-longitudinal resonance approximation |cos0| and the thermal cutoff for longi- 
tudinal (6 = 0°) propagation are also shown. 



754 



15° 2.022-* 

30° 2.074 - 

45° 2.136- 

60° 2.189 - 

75° 2.224 - 




FIGURE 4. The extraordinary mode index of refraction. 

The extraordinary-mode frequency band Fz < F < Fr- has no special designation nor reso- 
nance approximation. It is shown in figure 4 as a function of F with parametric and P = 2. Here, 
for a given F ^ P,\ cos | ^ cos Or. 

The polarization of these modes is quite similar. From (4.5) and (4.6) it can be shown that 
the electric vector of the whistler mode always rotates in the same sense as electrons gyrate about 
the static field. The rotation in the other mode has the same sense when F > P, but is opposite 
when F < P because ri- < 1. 

Although the indices of refraction given in (4.3) are complicated algebraically, when they are 
substituted into (3.23), the emission equation for both modes of propagation reduces to a quadratic 
in cos 2 0. v (F), 



^here 



A cos 4 8 + B cos 2 0* + C = O, 

A = j8|F 2 [(P2 - F 2 )» - F 2 (P* - F 2 )] + p 2 2 P 2 F 2 (F - s/y) 2 , 
B = /3 2 F 2 [2(P 2 - F 2 ) 2 + P 2 - 2F 2 ] (F - sly) 2 - P 2 (F - s/y)\ 

C = F\P 2 -F 2 + l)(F-sly)\ 



(4.11) 



(4.12) 



The allowable parameter ranges for /3 2 and P and allowable frequencies F are determined by the 
condition that the elementary roots, 



cos 2 Of = [-B±(B 2 - 4AC) l l 2 ]/2A, 



(4.13) 



are real, positive, ^ cos 2 Or according to the mode, and ^ 1. These superscripts which denote 
the two roots should not be confused with the subscripts (±) which denote the two modes of propaga- 
tion. In keeping with the emission directions assigned by (3.23), it is required that 0^0^^ 77-/2; 
^ 0f ^ tt/2 when F > sly, s > 0; or tt/2 ^ 6% ^ tt when F < 5/7, s > 0; and 0^6%^ tt/2 when 
5 < irregardless of F. 

In principle, one could determine F g (0), but the equation is an eighth degree polynomial and its 
analytical properties are more obscure. 

Since the two propagation modes are quite different, their emission properties will be discussed 
separately. 

Ordinary- (whistler-) mode emission is discussed first. A typical set of curves for cos 0^ (F) 
are plotted in figure 5 for s = 0,±l,±2,±3, with F = 2.25 and several arbitrary values of j8 2 . The 
Cerenkov emission is in the forward hemisphere, the normal cyclotron emission is in the backward 
hemisphere, and the anomalous cyclotron emission is in the forward hemisphere. The cos 
0^" and cos 0^ regions are adjacent but distinct with a boundary at cos 0'. defined below. In the 
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Figure 5. Solutions of the emission equation Of (F) for 
Cerenkov (s=0) and normal (s > 0) and anomalous 
(s < 0) cyclotron radiation in the ordinary (whistler) 
mode(P2/Pi = V3). 



Cerenkov case, only, there is an upper limit /3' 2 on (3 2 which occurs on the boundary cos 6' 8 ; however, 
there is always a root for small f$2 when the index is unbounded. It is apparent from the plots 
that a given frequency may be emitted along two cones defined by 0+ and 0~, along only one cone 
defined by 0f, or may not be allowed at all. An upper bound on the index which is discussed in 
the next section is also shown. 

The solutions of the emission equation for this mode have a few analytical properties which are 
quite easily demonstrated. In the low frequency limit F < 1 (P ~ 1, fc < 1), the solutions can be 
expanded in Maclaurin's series, 

cos 2 0:={(F-slyn(F~slyy-pi(2P^l)F 2 ] ipjF*(PlP* + s 2 ly 2 )} [1 - 0(/3|F 2 )], 

cos 2 6- = {cos 2 6 R I[1 - /3 2 (2P 2 + l)F 2 /(F- 5 /y) 2 ]} [1 + O^fF 2 )]. (4.14) 
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From these expressions it is easily verified that 1 > cos 2 6\ > cos 2 6 ~ > cos 2 Or, cos 2 0+ > 1 for 
5/0, and cos 2 d~ approaches cos 2 Or for 5 # 0. Hence, as F decreases, Fr+ decreases faster, so 
that n+ approaches the resonance limit at Fr+ and not the limit at F = 0. Also, the emission is 
transverse to B () when F <^ 1. In the Cerenkov limit = (5 = 0), the emission equation has the 
simple solution 



n(0)=|±|a- 



■WiPty 2 ) 1 ' 2 . 



(4.15) 



Since F must be real, the parameters must satisfy the condition fcPy ^ V2. In the cyclotron limit 
= or 7r, the emission equation reduces to an elementary cubic in F which always has one root 
F s (0 or 7r) in the range < F < 1. By definition, the boundary between the two solutions is 



cos 2 8+ = cos 2 0- = - B/2A = cos 2 0' , 



(4.16) 



and the vanishing discriminant B 2 = ^AC determines F' s and the maximum ($! 2 for the case 5 = 0. 

The frequency spectra and the angular patterns of the power given by (3.26) and (3.27) have 
been evaluated for a few specific parameters. The results for Cerenkov radiation and the first 
few cyclotron harmonics are presented in figures 6 and 7. The parameters of the background 
medium have been assigned the values fn = 0)bI2tt = 102 kc/s and fp = 0)pI2tt = 229 kc/s, which 
correspond to the magnetospheric locus L = 3, A = 30° defined in the appendix. The pitch angle 
1// between the electron's velocity and B has been set at 30° so that /WjSi = V3. In figures 6 and 7, 
the radiated power is shown for the electron kinetic energies KE= 10, 1000 keV (dashed and solid 
curves); the very fine structure involving AF < 0.01 is omitted. When both roots of the emission 
equation contribute to the power, the superscripts ± are explicitly given with the values of 5. The 
dotted curves are mean values of the power in regions with very large amplitude fluctuations where 
the index is unrealistically large. The effects of a bound on the index, which are discussed in the 
next section, are also included (solid curves only). 

These figures demonstrate several unanticipated features of the emission process. First, 
there is a broad maximum in the Cerenkov spectrum, and the maxima in the cyclotron spectra 
are not at Fi(tt) as linear emission theories tacitly assume. Second, the angular patterns show that 
the major portion of the energy is radiated at large emission angles from the static field B (wave 
normal direction, not ray path direction). Third, the energy in the cyclotron harmonics increases 
with harmonic number. Fourth, the energy radiated by the cyclotron harmonics becomes infinite 
as F— >0 or 0— »7r/2. Finally, the power radiated by a nonrelativistic electron increases as fit 
(or KE) decreases. (The infinite amount of energy indicated by these last three properties is sup- 
pressed by introducing thermal, ion, or collisional effects.) 
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FIGURE 6a. Frequency spectra of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s < 0) cyclotron radiation 
in the ordinary (whistler) mode from a 10 keV electron 
in a magnetoplasma. 

(1=3.0, \ = 30°) 




FIGURE 6b. Angular patterns of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s < 0) cyclotron radiation 
in the ordinary (whistler) mode from a 10 keV electron 
in a magnetoplasma. 

(L = 3.0, k = 30°) 
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FIGURE 7a. Frequency spectra of normal (s > 0) and 
anomalous (s < 0) cyclotron radiation in the ordinary 
(whistler) mode from a 1000 keV electron in a magneto- 
plasma. 

(L = 3.0, \ = 30°) 
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Figure 7b. Angular patterns of normal (s > 0) and anom- 
alous (s < 0) cyclotron radiation in the ordinary 
(whistler) mode from a 1000 keV electron in a magneto- 
plasma. 

(L = 3.0, A = 30°) 



These properties of the power expressions can be demonstrated approximately without per- 
forming extensive numerical analysis. From (4.5) and (4.6) it is easily verified that \ajc\ > \a$\ — 1 
for < 6 < rr. Using this result and the emission equation as a definition of n, and assuming that 
\Jg\ < \(XkJs\ and that the partial derivatives of n are slowly varying, the power expressions are 
proportional to 



dPoldf, d 2 P ldn ~ FVl tan 2 ^//3 2 3 (P 2 -F 2 ) 2 



(4.17) 



dPsldf, d 2 P s ldVi - s 2 FJ 2 s tan 2 0/y 2 /3 3 (P 2 - F 2 ) 2 s # 0, 



(4.18) 



for F # P. For small F the approximation cos 6~ — F verifies additional frequency properties. 
When 0=0 or 7T, it is found that a* = 0, ao = : +- 1, X= 0, and the power expressions vanish for all s 
except 



dPtldf= [ (7re 2 /c) (wffilfc) (1 + C F )I\1 ~ C e \ ]f-f x m. 
d 2 P 1 /dCL=[(e 2 l^7rc) (co 2 ^ 2 ) (1 + C f )I\1+PMI + Cf)\]f~f 1 m, 



(4.19) 



where 



C F = P 2 {2F - 1)/2(1 - F) (P 2 - F 2 + f), 

Ce = P 2 l2(l-F)(P 2 -F 2 ). 



(4.20) 



In addition to the results at L = 3, X = 30°, i// = 30°, the power expressions have also been 
evaluated at 3, 15°, 19°; 6, 45°, 19°; and 6, 15°, 6°, in order to ascertain the effect of varying /b,/p, 
and i//. Although the relevant frequencies vary considerably, the general shape of the frequency 
spectra and angular patterns do not change appreciably. As a quantitative example, it is found 
that a 10 keV electron at 6, 15°, 6° radiates less than 10~ 33 W/(c/s) at frequencies above 1 kc/s (the 
lowest considered) although/ fi = 5.5 kc/s and/p = 70 kc/s so that n is large, and furthermore, below 
1 kc/s the cyclotron emission is confined to 90° < < 100°. 

Extraordinary-mode emission is considerably more complicated than ordinary mode emission. 
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FIGURE 8. Solutions of the emission equation Of (F) for 
Cerenkov (s = 0) and normal (s > 0) and anomalous 
(s < 0) cyclotron radiation in the extraordinary mode 
(j8 2 /j3i = V3). 
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The emission equation roots cos Of{F) are plotted in figure 8 for 5 = 0, ±1, ±2, ±3, ±4, with 
P = 2.2S and a few arbitrary values of fa- The Cerenkov emission is again in the forward hem- 
isphere, but the cyclotron emission is spread over both the forward and backward directions. The 
complicated nature of the solution for F ~ sly ^ P is qualitatively expected from (3.23), and the 
discontinuity at F = P for small fa is attributable to the sharp change in n~ there. 



Emission in this mode is restricted by the index to the frequency band F/<F<F^ AX 
= (P 2 +iy 12 . Some analytical properties of the solutions are observed by expanding them in 
the limit fa « 1, 



cos 2 0+ = (1 - slyFfiPl + 8 + (P 2 - F 2 )/P 2 - (P 2 - F*flP*F* + 0(/3 2 ), 
cos 2 0- = cos 2 ^[l-)3 2 8/(l-5/yP) 2 + 0(/3|)], 



(4.21) 
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fhere 



8 = 2F 2 IP 2 - 1 - 2(P 2 - F 2 ) 2 /P 2 - cos 2 6 R . 



(4.22) 



In general, F =f s/y so that cos 2 0+ > 1 due to the first term in its expansion and cos 2 B~ — cos 2 
6r{< 1 for F > P). When F — s/y, the root cos 2 6 S + may contribute to the emission pattern depend- 
ing on the other parameters. However, the expansion for cos 2 0~ is not valid when F — s/y. At 
F = P 7^ s/y, it is found that 8 = and cos 2 6r—1 so that cos 2 6~ — 1. Also, at the frequency band 
limit F = (P 2 + 1) 1/2 , cos 2 8' - cos 2 Or = 0. 

The frequency spectra and angular patterns of the Cerenkov and cyclotron power radiated 
in the extraordinary mode by 10 and 1000 keV electrons at L = 3, A. = 30°, with i// = 30° are shown 
in figures 9 and 10. In this mode much of the energy is radiated at small angles to the static field 
Bo by the cyclotron fundamental. The power singularity at # = 7r/2, F = (P 2 + 1) 1/2 is suppressed 
by a bound on the index as before. The Cerenkov emission near = 0, F = P is not adequately 
described by this formalism; it consists of electric waves (E x , Ee — 0, Ek^O) which travel with 
the source (J3 2 n = 1) and the infinite power obtained here is invalid. A detailed treatment of this 
Cerenkov effect has been given by Gerwin and Guernsey [1963]. 
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FIGURE 9a. Frequency spectra of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s > 0) cyclotron radia- 
tion in the extraordinary mode from a 10 keV electron 
in a magnetoplasma. 

(L - 3.0, k = 30°) 



FIGURE 9b. Angular patterns of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s < 0) cyclotron radia- 
tion in the extraordinary mode from a 10 keV electron 
in a magnetoplasma. 

(L = 3.0, A = 30°) 
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FIGURE 10a. Frequency spectra of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s < 0) cyclotron radia- 
tion in the extraordinary mode from a 1000 keV electron 
in a magnetoplasma. 

(L = 3.0, a. = 30°) 



FIGURE 10b. Angular patterns of Cerenkov (s = 0) and 
normal (s > 0) and anomalous (s < 0) cyclotron radia- 
tion in the extraordinary mode from a 1000 keV electron 
in a magnetoplasma. 

(L = 3.0, A. = 30°) 
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5. Discussion 

Since the elementary model plasma described in the last section leads to the emission of an 
infinite amount of energy by the gyrating electron, it is necessary to introduce some additional 
effects which will limit the emission and give a physically meaningful result. The major factors 
which were omitted previously are thermal motion, ionic motion, and collisions in the background 
medium. For VLF propagation in the magnetosphere and certain other applications, the dominant 
modification of the propagation properties is obtained by including thermal motion alone. 

Mathematically, the boundless energy radiated into the plasma medium is due to the un- 
bounded property of n± at Fr±. From figures 5 and 8 it is evident that the solution of the emission 
equation approaches Or when F approaches 0, approaches 7r/2, or \s\ increases, and the power 
increases sharply because otk is proportional to n 2 . Consequently, a physical process which 
introduces a finite upper bound n ± on the index is essential if the radiation is to have physical 
meaning. Furthermore, such a bound implies |cos 0~|^|cos 0±|, F^F±, |s|<F±, and fi>> > fiz+ 
according to (3.23), and the values of these limits depend on the explicit form of n±. 

An artificial bound on the index is introduced here to demonstrate the qualitative effect that 
it produces on the frequency spectra and angular patterns of the radiated power. The index is 
assumed to be given by (4.3) when n ± < n ± , and by n± when n ± > n ± . Since n± increases mo- 
notonically as |cos 6~\ approaches cos Or, the solution of the emission equation is cutoff at the 
discontinuity where |cos #~|=|cos 6±\. In figures 5 to 10 the cutoff is shown for the case n = 20. 
The finite power for this artificial case is described by the solid curves. Although the sharp 
edges would disappear for a continuous index, the dominant maxima are expected to remain. It is 
noteworthy that the nonrelativistic maxima in both propagation modes occur at the cutoff (figs. 6 
and 9) so that a small change in the cutoff location produces a large change in the radiated power. 
The extent to which this artificial case is quantitatively applicable to a real case depends on the 
validity of the index approximation which is discussed below. 

Insofar as the index approximation is valid, a comparison can be made between the power 
in the different harmonics and between different electron energies. For this purpose the power 
expressions were evaluated at a few intermediate energies (not shown). In general, the maximum 
power in the normal and anomalous cyclotron harmonics, and the Cerenkov radiation (when it 
exists) is comparable at a given energy. The total radiated power is a slowly decreasing function 
of fit (or kinetic energy for a fixed pitch angle) at nonrelativistic speeds, with an average level of 
10 -30 W/(c/s) in the ordinary mode and 10~ 25 W/(c/s) in the extraordinary mode. The range of 
speeds y8 2 , which allow Cerenkov radiation, is severely restricted from below by a bounded index 
and above by lack of a polarization interaction. At relativistic speeds the higher cyclotron har- 
monics dominate, but they are restricted by the upper bound s". 

The formulas for the dielectric tensors, polarization coefficients, and refractive indices which 
include additional effects are quite complicated and difficult to apply. The extent of their influence 
on the results is evaluated approximately here by considering some special cases. 

Random thermal motion is considered first. Since it tends to disorganize the ordered motion 
which is essential for wave propagation, there is a minimum phase velocity (or n ± ) which the 
thermal medium can support. The dielectric tensor for a thermal medium is known, and an ex- 
tensive review of the theory has been given by Stix [1962]. However, the dispersion equation 
is so complicated, in general, that it has only been solved in a few special cases. In particular 
for 0~0 (or 7r), the maximum occurs at F = l, and for a Maxwellian distribution with temperature 
T, its value is K + - ?> 1 ' 2 P 2 ^{ttUc 2 IKT)^I2. Liemohn and Scarf [1963] used a bell-shaped alge- 
braic distribution which gives a much simpler dispersion equation to study the longitudinal whistler 
mode n+(0 = 0). They found that the zero-temperature index is a valid approximation when 
1 — F » {KT/Mc 2 y l2 n+. As yet n+ has not been studied for arbitrary 0, and n- has not been con- 
sidered at all, although similar behavior is anticipated. 

It is impossible to accurately assess the thermal effect on the emission with this limited 
amount of information about the index. Nevertheless some speculation is conceivable. For the 
magnetospheric case T= 1250 °K and P = 2.25, the zero-temperature index for = or tt is valid 
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when 1— F » 0.02 and n + ^ n + = 20 (see fig. 3). Hence, the artificial discontinuous index used 
in figures 5 to 7 is a good approximation when |cos 0\ >> 0. At other angles of emission and in 
the extraordinary mode the validity of the index approximation is unknown at present. 

In a thermal magnetoplasma the electromagnetic waves are also absorbed as they propagate so 
that the power expressions are only valid near the source. The absorption is due to the cyclotron 
resonance mechanism. The component of thermal motion parallel to B causes a Doppler shift 
of the propagation frequency in the rest frame of each electron. For each angle of propagation 
there is a corresponding velocity component which shifts the propagation frequency to the reso- 
nance. Hence, the damping decrement depends on the frequency and angle of propagation as 
well as the thermal velocity distribution of the electrons. Its value is determined by the general 
solution of the dispersion equation. Due to the frequency and angular dependence, the power 
spectra and patterns will change as the signal propagates away from the source. Fortunately, 
in many applications the temperature of the plasma is so low that very few electrons resonate and 
the thermal attenuation is negligible. If the background plasma has an appreciable number of 
"fast" electrons due to a high temperature or other sources, however, the damping must be 
considered. 

In a recent series of papers, several limiting cases of emission in a thermal magnetoplasma 
are analyzed by Pakhomov, Stepanov, and Aleksin [1962], Pakhomov and Stepanov [1963a, 1963b]. 
In general, their analysis is restricted to the index values n± ~ 1 where several approximations 
for the emission equation and the power could be made. When the limitation on n± is relaxed, 
other restrictions on o> or 6 are imposed. Hence, their results cover a variety of special cases 
but do not provide information about the complete frequency spectra and emission patterns. 
A convenient feature of their expressions is that the thermal attenuation factors are included 
explicitly. 

Due to the inertia of the ions, they only interact appreciably with the electromagnetic fields 
when F « 1. In this region Hines [1957] has shown that the index is bounded for 6 ~ 6r and 
ionic effects are important when F < (m/M) 1/2 . Since cos 6~ ~ F when P > 1 according to (4.14), 
the limit is also applicable to the angle of emission. For the particular case of ionized hydrogen, 
the limit is so close to the transverse direction that the bound imposed on the power is extremely 
large compared to the anticipated thermal bound. Consequently, ionic effects can usually be 
disregarded. 

Collisions which the electrons make with heavy particles in the medium are usually approxi- 
mated by a viscous damping force in the equation of motion. Such a force restricts the wave 
propagation in a manner similar to thermal motion. The theory has been reviewed by Ratcliffe 
[1959], Ginzburg [1961] and others. In general, collisions are only important when the collision 
frequency o) c satisfies o) c > o). In the magnetosphere, o) c < 10 rad/sec so that collisions are neg- 
ligible in the VLF band. Unfortunately ,the dielectric tensor for this elementary theory is not 
hermitian so that when collisions must be included a different formulation of the problem is needed. 

6. VLF and LF Emissions 

The foregoing properties of the radiation from electrons in a magnetoplasma suggest some 
conclusions about the role that electron emissions might play in the generation of VLF and LF 
emissions. An important test for any prospective source is to compare the theoretical and ob- 
served power of the signal. Also, since the emissions have a wide variety of discrete dispersion 
patterns as well as continuous bands, some mechanism is needed to provide for frequency selec- 
tion. Finally, the emission properties suggest another strong signal which is only observable by 
satellites of rockets. 

In the magnetosphere, the geomagnetic field and the plasma density vary so slowly in space 
and time that the propagation theory for uniform media is usually acceptable. However, the 
magnetosphere changes significantly over several wavelengths so that a signal which propagates 
in one region may be absorbed at a resonance or reflected at a cutoff in a nearby region. A further 
complication is presented by the ionosphere where the electron density increases sharply. On the 
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geomagnetic equator the cyclotron frequency varies continuously as the inverse cube of geo- 
centric distance from ~ 1 Mc/s at the earth's surface to ~ 1 kc/s near 10 earth radii. The variation 
of the plasma frequency is very similar with a range from ~ 5 Mc/s in the F2 peak of the ionosphere 
to —20 kc/s near 10 earth radii. Mathematical models of the region which include latitude 
effects are presented in the appendix. In the magnetosphere, P varies from ~ 0.03 to ~ 20. 

At the conclusion of section 3, it was noted that the ray path followed by the electromagnetic 
energy is not in the wave normal direction (#), but rather is in the time-average Poynting vector 
direction much closer to the magnetic field direction. As a result, it is frequently assumed that 
the ray path follows a magnetic flux tube. While such an approximation is questionable, it provides 
a convenient framework for discussing some general propagation properties of electron emissions 
from their sources to potential receivers. 

Due to the cutoffs and resonances of the modes, the signals only propagate in well defined 
regions of the magnetosphere. In the whistler mode the emitted frequencies can always propagate 
into regions of increasing cob (inward), but can only propagate into regions of decreasing cob (out- 
ward) if F < Fr+ (neglecting thermal and collisional effects). In the extraordinary mode the 
emitted frequencies can propagate into regions of increasing a># until F = Fz, where they are 
reflected and can propagate into regions of decreasing lob if F < Fr— Hence, the signals in the 
whistler mode can penetrate the ionosphere and reach the ground, whereas signals in the other 
mode are reflected, in general, above the ionosphere, and both modes have a minimum upper 
frequency limit which occurs at the equator. A schematic diagram of possible propagation paths 
is shown in figure 11. Since part of the whistler mode energy is reflected at the ionosphere by 
irregularities, an emission may echo between hemispheres. Similar echoing in the extraordinary 
mode is only possible when the source is very near the equator. 

A detailed calculation of the VLF and LF power which is radiated into the magnetosphere is 
impossible at this time because the propagation properties are not fully understood and the phase 
space distribution of the electrons is unknown to a large extent. Hence, some crude estimates of 
the power from incoherent and coherent sources must suffice. 

First, it will be shown with a specific numerical calculation that incoherent electron emissions 
apparently cannot account for the observed emissions received on the ground in the whistler mode. 
The density of magnetospheric electrons N ( > which emit appreciable radiation is assumed to be 
~ 0.1 elect/cm :j , which is fairly consistent with satellite measurements reported by O'Brien [1962] 
and whistler analysis by Liemohn and Scarf [1964]. The volume of a flux tube with a 1 cm 2 
base which crosses the equator at 3 earth radii (L = 3) is found to be ~ 10 10 cm*'*. On the basis 
of the numerical results for the ordinary mode, the average emitted power P e along the flux tube is 
assumed to be on the order of ~ 10~ :l ° W/(c/s). With these gross approximations the incoherent 
power at the base of the flux tube is found to be 10 -21 W/cm 2 (c/s). This is several orders of mag- 
nitude below the observed power level P r of 10~ 14 W/cm 2 (c/s) at the base of the magnetosphere 
(above the ionosphere). 
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Figure 11a. Possible ray (energy) paths for the ordinary 
(whistler) mode of propagation in the magnetosphere. 



Figure lib. Possible ray (energy) paths for the 
extraordinary mode of propagation in the magnetosphere. 
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Second, coherent emissions from "bunches" of electrons are shown to be a feasible source 
of the emissions. Since the electromagnetic fields of individual sources must be in phase for 
coherent radiation, the total power is proportional to the number of emitters squared. The 
required size of the bunch is determined by equating the observed power over area A r to the 
emitted power from an extended source of cross section A e and length l e with coherent efficiency e. 

A r Pr = Pe(eN e Aele) 2 . (6.1) 

The radiation is assumed to be confined to a flux tube so that B r A r = B e A e . Then, with the values 
of P r , P e , and N e given above and the arbitrary assumptions B e IB,=0.1 and € = 0.001, eq (6.1) 
reduces to A e le = 10 +23 cm 4 . Hence, the dimensions of the bunches are small relative to the 
magnetosphere, and the scale is comparable to a few wavelengths. On the basis of this crude 
calculation, it appears feasible that coherent emission from electron bunches is a conceivable 
source of the emissions in the ordinary mode. 

Since the electron emissions are continuous over a wide band of frequencies, the discrete 
nature of the observed VLF emissions must be attributed to the coherent interaction. The trig- 
gering of emissions, which Helliwell [1963] has observed, suggests that strong whistler mode signals 
from lightning strokes interact with the electrons and organize their phase space distribution so 
that specific frequencies are amplified by coherence. The various interaction mechanisms which 
may produce this bunching have been reviewed by Brice [1964] . 

The one-dimensional calculations by Dowden [1962c] and others, which support the theory 
that emission at the Doppler-shifted fundamental cyclotron frequency (^1(77)) is responsible for 
hooks, is not ruled out by these results, but it is much more likely that stronger emissions at other 
frequencies in nonlongitudinal directions dominate the coherent signals. The traveling wave 
tube amplification of incoherent signals which is proposed by Gallet and Helliwell [1959] seems 
improbable because a gain of several orders of magnitude is required. 

The power in the extraordinary mode is considerably stronger than that in the ordinary mode 
above the ionosphere. Hence, much less coherence is necessary for this mode to produce signals 
which are detectable in the presence of the observed whistler mode noise level. The average 
power level is on the order of 10 -25 W/(c/s) so that 10 9 electrons in the flux tube described above 
produce 10 -16 W/cm 2 (c/s) of incoherent power at the base of the tube. Rocket observations by 
Walsh, Haddock, and Schulte [1963] have revealed intense noise signals in the extraordinary 
mode so that this mode cannot be ignored. 

Certainly the calculations presented here are based on very limited numerical results and must 
be considered as only preliminary estimates of the power levels. A thorough quantitative analysis 
must await the inclusion of additional limiting physical effects such as thermal motion and perhaps 
collisions. However, they serve to point out some inherent weaknesses in the existing theories 
for VLF emissions in the ordinary mode, and they point out complications by the extraordinary 
mode when observations are made above the ionosphere. Furthermore, order of magnitude cal- 
culations of the power provide an additional test for future theories of VLF and LF emissions which 
are based on charged particle radiation in a plasma. 
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7. Appendix 

The magnetic field of the earth is closely approximated by a dipole field of the form 

B = 6(-2 sin AJJ + cos k 1)1 R\ 

where R is the geocentric radial distance and k is geomagnetic latitude. With R expressed in 
earth radii the dipole moment b is 0.313 gauss (earth radii) 3 . The field line equation is 

R —L cos 2 \, 

where L is the radial distance at which the field line intercepts the equator (Mcllwain L parameter). 
The scalar field is given by 

£o = 6(l + 3 sin 2 k)^IU cos 6 k. 

Several mathematical models which closely approximate the thermal plasma density in the 
magnetosphere have been proposed. In this analysis the field-line diffusion model by F. S. Johnson 
[19621 is used: 

/Vo= 12500/?-* cos 6 k exp (GR l sin 2 k) R > 1.3 

where G = MgR r j2KT is the scale factor of the medium in earth radii (R ( ). The temperature of the 
plasma Tis usually assumed to be ~ 1250 °K, which gives G = 3. 

The motion of high-energy electrons in this medium is described by the first-order orbit theory 
for trapped particles [Welch and Whitaker, 1959; Chamberlain, 1963]. The electrons are con- 
strained to follow helical orbits along magnetic flux tubes and are reflected at conjugate mirror 
points in opposite hemispheres. The longitudinal motion is neglected here. Due to conserva- 
tion of the electron's magnetic moment, the pitch angle \\f between the velocity and field vectors 
is related to its equatorial value \fj e by 

sin t// = (l + 3 sin 2 k) l/4 sin i/^/cos 3 k. 



The latitude of the mirror point k m is determined by setting i// = 7r/2 in this equation. Hence, the 
orbit is completely specified by only two parameters, L and i// r , which are independent of the kinetic 
energy. 
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